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Abstract. New finite energy sum rules (FESR's) for extracting m s from hadronic r 
decay data are constructed which (1) significantly reduce potential theoretical uncer- 
tainties present in existing sum rule analyses and (2) remove problems associated with 
both the poor convergence of the OPE representation of the longitudinal part of the 
us vector and axial vector correlators and the large statistical errors in the us spectral 
data above the K* region. 



The ratio of the hadronic r decay rate through the / = ij — ud, us vector (V) 
or axial vector (A) current to the corresponding electronic decay rate, R^/A'M = 
r[r~ — > v T hadronsy/A;y (7)]/r[ r_ — > u T e~u e ( , y)), can be written [1]: 
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with y = s/m 2 , Vij the / = ij CKM matrix element, Sew & n electroweak correction, 
and PyJ A .ij{s) the spectral function of ^-v/A;ij( s ) > ^v/Aij( s ) (J = 0,1) being the 
spin J part of the / = ij vector (V) or axial vector (A) correlator. The second 
line follows from the first as a consequence of the general FESR relation, valid for 
any II without kinematic singularities, and any w(s) analytic in the region of the 
contour, f*° ds w(s)p(s) = ^ §\ s \=s dsw(s)H(s). Experimental data thus allows 
access to the ud-us spectral difference, and hence to integrals of the corresponding 
correlator difference, which, for large enough sq, are dominated by the D = 2 term 
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of the OPE representation, proportional to m 2 . For such s one may, therefore, 
hope to extract m s in terms of appropriately weighted integrals of the experimental 
spectral data. 

A complicating factor in any attempt to determine m s based on this observation 
is the non-convergence of the OPE representation of the longitudinal ((J) = (0)) 
integral at scales < m 2 [2]. The current inability to make an experimental longitu- 
dinal/transverse separation above s ~ 1 GeV 2 thus, at present, precludes a reliable 
analysis using sum rules with significant longitudinal contributions. Recent analy- 
ses [3,4], which either work with the spectral data without making a longitudinal 
subtraction, or attempt to place loose experimental bounds on the longitudinal 
contribution, employ "spectral weights" (for the transverse ((0 + 1)) case, defining 
y = s/sq, these are w*(y) = [(1 — y) 2 (l + %)](1 — y) k )- 111 this paper we work with 
the "transverse" (V + A)^ 1 ^ — (V + ^4)„° s +1 - ) difference, and try to find alternate 
weight choices which improve the reliability of the analysis. The D = 2 term in 
the corresponding OPE is 



where a = a s (Q 2 )/ii. Details of the treatment of D = 4, 6, 8 contributions may be 
found in Ref. [5]. 

On the spectral side we employ the ALEPH ud and us distributions [3,6]. The 
longitudinal subtraction is performed using sum rule methods [5]; for the weights 
employed in our analysis the integrated longitudinal subtraction represents < 1% of 
the integrated D = 2 OPE ud-us difference, making the impact of any uncertainties 
associated with this procedure negligible. 

The weights employed in the ud-us FESR's of our analysis have been chosen so as 
to reduce both theoretical and experimental difficulties. On the experimental side, 
we seek to (1) de-emphasize us spectral contributions from the region above the 
K*, since the ALEPH determination of pv+A-,us has ~ 20 — 30% statistical errors 
in this region [3] and (2) reduce the strong cancellation in the ud-us difference, 
which otherwise greatly magnifies the impact of experimental uncertainties. (See 
Ref. [5] for a detailed discussion of the second point.) Weights which fall more 
strongly with s above s ~ 1 GeV 2 decrease the level of ud-us cancellation and 
simultaneously suppress high-s contributions, decreasing the impact of both the 
longitudinal subtraction and experimental errors on pv+A-,us- 

On the theoretical side, the goal is to control an important potential theoretical 
systematic uncertainty. For the weights the contour improvement prescrip- 
tion [7] is known to produce a significant improvement in the convergence of the 
known terms of the integrated D = 2 series. The smallness of the last (0(a 2 )) 
known term, however, turns out to result from strong cancellations between contri- 
butions from different regions of the circular part of the contour, \s\ = s [5]. Since, 
assuming continued geometric growth of the coefficients, similar cancellations do 
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FIGURE 1. The agreement between the OPE and hadronic sides of the FESR corresponding 
to the weight, w 2 o(y) for 2.55 GcV 2 < s < m\. The solid line is the OPE side, using the values 
of m s and the D — 8 contribution obtained in the fitting procedure described in the text and, in 
more detail, in Ref. [5]. The dashed line is the hadronic side, obtained using the ALEPH spectral 
data from which the longitudinal component has been subtracted as described in Ref. [5]. 



not persist to higher orders [5], an estimate of the truncation error based on the 
size of the 0{a 2 s ) term is unreliable. We have constructed 3 alternate polynomial 
weights which both avoid such "accidental" cancellations and emphasize regions of 
the plane for which the convergence of the D = 2 series is optimal. The result is 
a very strong suppression of possible higher order D = 2 contributions [5]. The 
explicit forms of the weights, as well as details of this improvement, are given in 
Ref. [5]. 

The results of our analysis are as follows. First, all 3 new weights yield consistent, 
and stable, values of m s in the window 2.55 GeV 2 < so < m r- An illustration of this 
fact is given, in Table 1, for the weight tuio(y) = [1 — 2/] 4 [l + 2/] 2 [l + 2/ 2 ][l +U + y 2 ] = 
1 — y — y 2 + 2y 5 — y 8 — y 9 + y 10 , where y = s/sq, which is favorable from a theoretical 
point of view because the absence of y 3 , y A terms removes D — 8, 10 contributions 
to the integrated OPE. Second, for our two other weights, which do not share this 
property, the D = 8 contributions, which are determined self-consistently, are also 
stable in this window. These conditions are not satisfied for the FESR's based 



TABLE 1. The extracted value of m s (l GcV 2 ) in MeV as a 
function of s for the weight Wio having no D = 8, 10 contri- 
butions. 

s (GeV 2 ): 2.35 2.55 2.75 2.95 3.15 

m s \l GeV 2 ) (MeV): 153.2 159.0 162.2 163.4 163.2 

on the spectral weights, w*. We choose, for our final analysis, that weight among 
the three constructed above (called u>2o in Ref. [5]) which leads to the smallest 
fractional statistical error. The match between the OPE and hadronic sides of the 
corresponding FESR which results once m s and the D = 8 contribution have been 
optimized is shown in Figure 1, and is clearly excellent. Our final numerical result 
for m s , in the MS scheme, is 

771,(1 GeV 2 ) = 158.6 ± 18.7 ± 16.3 ± 13.3 MeV , (3) 

which, using four-loop running, corresponds to 

m s (A GeV 2 ) = 115.1 ± 13.6 ± 11.8 ± 9.7 MeV . (4) 

The first error is statistical, the second due to the uncertainty in \V US \, and the 
third theoretical, with the latter dominated by our estimate of the error associated 
with truncating the D = 2 series at 0{a 2 s ). Improvements in the accuracy of the us 
spectral data, such as will be possible at BaBar, will serve to significantly reduce the 
first error and, simultaneously, allow use of weights other than w 2 o which produce 
reduced theoretical truncation errors. 
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